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We report measurements of the bulk magnetic susceptibility and 29 Si nuclear magnetic resonance 
(NMR) linewidth in the heavy- fermion alloy CeRhRuSi2- The linewidth increases rapidly with 
decreasing temperature and reaches large values at low temperatures, which strongly suggests the 
wide distributions of local susceptibilities \i obtained in disorder-driven theories of non-Fermi-liquid 
(NFL) behavior. The NMR linewidths agree well with distribution functions P(x) which fit bulk 
susceptibility and specific heat data. The apparent return to Fermi-liquid behavior observed below 
1 K is manifested in the vanishing of P(x) as X ~~ * 00 1 suggesting the absence of strong magnetic 
response at low energies. Our results indicate the need for an extension of some current theories of 
disorder-driven NFL behavior in order to incorporate this low-temperature crossover. 

PACS numbers: 71.27.+a, 75.30.Mb, 76.60.Cq. 



I. INTRODUCTION 

A breakdown of the standard Landau Fermi-liquid the- 
ory is signaled in certain heavy-fermion metals by anoma- 
lies in thermodynamic, transport, anoLoptical properties 
at low temperatures and frequencies a Although excep- 
tions exist, the anomalous properties are usually as fol- 
lows: the Sommerfeld specific heat coefficient j(T) — 
C(T)/T diverges as — lnT; the magnetic susceptibil- 
ity x(T) varies as 1 — aT 1 / 2 or diverges as — lnT or a 
weak inverse power of temperature; the electrical resis- 
tivity departs linearly with temperature from its T = 
value; and optical conductivity experiments in the non- 
Fermi-liquid (NFL) alloy UCu3.5Pdi.5 indicate a trans- 
port relaxation rate which varies linearly with frequency 
at low temperaturesB 

Attempts to understand this NFL behavior invoke one 
or more characteristics common to most such systems, 
viz., the possibility of an unconventional Konda effect pel 
proximity to ^quantum critical point (QCP),ETlI struc- 
tural disorder Jl or a combination of the latter twoB Re- 
cent experimental work has stressed the role that disor- 
der can play. In particular, the observed inhomogeneous 
broadening of copper nuclear magnetic resonance (NMR) 
lines in the NFL alloys UCus-^Pdz, x = 1.0 and 1.5, 
could be described by a disordep-iiiduced spatial distribu- 
tion of local susceptibilities Xj &3x3 Such a susceptibility 
distribution originates in the interplay between structural 
disorder and many-body effects intrinsic to /-electron 
systems, such as the Kondo effect and the Ruderman- 
Kittel-Kasuya-Yosida (RKKY) interaction between the 
magnetic moments. 



It is clear that if an interaction between moments is 
present the term "non-Fermi liquid" must be used with 
care, since strictly speaking Fermi-liquid theory deals 
only with the lowest-lying excitations of a system of in- 
teracting fermions and hence is correct only in the zero- 
temperature limit in the absence of a phase transition. 
Any kind of magnetic phase transition or glassy spin 
freezing at nonzero temperature invalidates this condi- 
tion, and NFL behavior is no longer a surprise. Following 
convention in this field, we nevertheless continue to des- 
ignate as NFL systems materials in which the properties 
mentioned above are found at intermediate temperatures, 
with the proviso "nearly NFL" applied if there is evidence 
for a crossover to a new state at low temperatures. 

In this paper we report measurements of the magnetic 
susceptibility and 29 Si NMR linewidth in the nearly- 
NFLE 2 ] heavy-fermion alloy CeRhRuSi2 , and consider the 
data in the light of two such disorder-driven scenarios: 

(1) _the so-called "Kondo disorder" picture ofJBernal et 
al.tB and Miranda, Dobrosavljevic, and Kotliar,B in which 
the RKKY interaction is disregarded and the local sus- 
ceptibility distribution is associated with a corresponding 
distribution of single- ion Kondo temperatures Z*-, and 

(2) a recent model by Castro Neto and co- workers! -based 
on the existence of quantum Griffiths singularities^ in a 
disordered-jsystem with RKKY couplings which is close 
to a QCPBi In the latter case Kondo and RKKY phe- 
nomena compete with each other in the random environ- 
ment, and the susceptibility is associated with fluctua- 
tions of magnetic clusters!] Both of these models have 
been shown to account for the obsefvwd susceptibility 
and NMR broadening in UCu 5 _ a Pd x Ja'E] and the Kondo- 



disorder model is in agreement with the transport and op- 
tical data for \]Cu5- x Pd x alloys. The Griffiths-phase pic- 
ture describes the thermodynamic properties of a number 
of NFL materials^ 

It should be pointed out that N-EL mechanisms based 
on unconventional Kondo effectsErQ have to date been 
treated only for isolated / ions. It would be useful to 
extend such pictures, first to the case of ordered /-ion- 
based compounds and then with the inclusion of disorder 
in analogy with the Griffiths-phase model. 

Similarities and differences between the Kondo-dis- 
order and Griffiths-phase theories of NFL behavior and 
the necessity for their modification at low temperatures 
are discussed in the light of our experimental data. Fits of 
the theories to the temperature dependence of the bulk 
susceptibility determine the parameters of each model, 
each of which then predicts the temperature dependence 
and size of the NMR linewidth with no further adjustable 
parameters. The measured linewidths are in good agree- 
ment with both. models. This corroborates the conclusion 
of Graf et al.]El based on susceptibility and specific heat 
measurements, that disorder-driven NFL behavior is im- 
portant in this system. 

The isostructural alloy system Ce^ui-z Rh a )2Si2 ex- 
hibits a variety of complex behavior associated with the 
Kondo effect and magnetio4pteractions.Ej t3 The phase 
diagram of this systemOlia is shown in Fig. [I]. The 
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FIG. 1. Magnetic phase diagram of Ce(Rui_ :E Rh a ;)2Si2. 
Data from Refs. [16|-]l8[ Paramagnetic (PM) and antiferro- 
magnetic (AF) regions are indicated. Circles: AF ordering 
temperatures Tjv. Triangles: Kondo temperatures Tk- The 
curves are guides for the eye. The arrow indicates the con- 
centration for which NFL behavior was observed by Graf et 
al. (Ref. 

end compound CeRu2Si2 is a heavy-fermion metal which 
shows no long-range magnetic order dawn to 20 mK,E£l 
whereas CeRli2Si2 is an antiferromagnetO with Neel tem- 
perature Tjv = 36 K. Antiferromagnetism is found for 
low (0.1 < x < 0.3) and high (0.5 < af-S^ 1) rhodium dop- 
ing. Neutron diffraction experimentst3 in the low-doping 



range show that the antiferromagnetism is incommensu- 
rate (i.e., the 4/ electrons are itinerant), but becomes 
commensurate, indicative of local moments, for x > 0.5. 
For x < 0.4 the specific heat indicates a characteristic 
energy scale, usually associated with the average Kondo 
temperature Tk of the material, which is higher than Tn 
in this composition range. 

The concentration x = 0.5 is near the critical value for 
suppression of to zero. In spite of its stoichiomet- 
ric composition CeRhRuSi2 is a disordered alloy, as the 
Rh and Ru atoms occupy the same crystallographic site 
and there is no evidence of superlattice formation. For 
x = 0.5 Graf et al£3 found the weak divergences charac- 
teristic of NFL behaviora in j(T) and x(T) for temper- 
atures between 1 and 30 K. There was no evidence for 
magnetic ordering, and Graf et al. concluded that NFL 
phenomena in CeRhRuSi2 are driven by structural disor- 
der. To our knowledge the region 0.3 < x < 0.5 has not 
been examined for NFL behavior. 

Below 1 K j(T) was seen to saturate, suggesting that 
CeRhRuSi2 exhibits a crossover from a region of anoma- 
lous magnetic response to a Fermi-liquid ground state as 
the temperature is lowered. It should be noted, however, 
that recent specific heat and ac susceptibility studies of 
UCu5_ K Pd x EJcj suggest that saturation of j(T) in this 
system may be associated with magnetic ordering, pos- 
sibly of a spin-glass nature or in the form of superpara- 
magnetic clusters. More information on the behavior of 
NFL systems at low temperatures is clearly needed. 

The bulk susceptibility agrees well with the Kondo- 
disorder model but is overestimated by the Griffiths- 
phase picture at low temperatures. This is not surprising, 
since the possibility of nearly-NFL behavior is built into 
the Kondo-disorder model, whereas the Griffiths-phase 
theory in its present form neglects effects, such as trans- 
verse Kondo fluctuations or residual interactions between 
clusters, which could modify or remove the Griffiths sin- 
gularities which cause the NFL behavior. 

The experimental data obtained to date do not dis- 
criminate clearly between the two pictures, since the dis- 
tribution function P{x) which describes the inhomogc- 
neous distribution of susceptibilities contains no qualita- 
tive feature sensitive to the existence of RKKY-coupled 
clusters. We speculate that dynamical properties such 
as nuclear spin-lattice relaxation rates may be more sen- 
sitive to low-lying excitations associated with spin-spin 
couplings, particularly at low temperatures, and suggest 
that further measurements of spin relaxation rates be 
carried out. 

bee. H of the paper describes our measurements of 
bulk magnetic susceptibility and 29 Si NMR spectra in 
CeRhRuSi2- The relation between inhomogeneity in 
the susceptibility and the NMR linewidth is reviewed in 
Sec. |ffl|. Sec. treats the single-ion Kondo-disorder and 
Griffiths-phase disorder-driven NFL mechanisms. Anal- 
ysis of the susceptibility and NMR data is discussed in 
Sec. and Sec. [vl gives our conclusions. 
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II. EXPERIMENT 

Measurements of bulk susceptibility and 29 Si NMR 
spectra were carried out on unaligned and field-alignecO 
powder samples of CeRhRuSi2 at a frequency of 20.220 
MHz and temperatures in the range 4.2-230 K. Field- 
swept NMR spectra were obtained using pulsed-NMR 
spin-echo signals and the frequency-shifted-and-summcd 
Fourier-transform processing technique described by 
Clark et al£3 The solid curve in Fig. || shows a 29 Si 
NMR spectrum from an unaligned powder sample at 4.2 
K. We attempted to fit this spectrum to an anisotropic 
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FIG. 2. Field-swept 29 Si NMR spectra in CeRhRuSi 2 for 
T — 4.2 K and spectrometer frequency 20.220 MHz. Solid 
curve: unaligned powder sample. Dashed curve: field-aligned 
powder, Ho || c. Dash-dot curve: field-aligned powder, 
Hole. 

powder patternE! convoluted with a Gaussian broadening 
function, but found a poor fit if the width of the broad- 
ening is assumed independent of crystallite orientation. 
The low-field side of the spectrum, which is due to those 
crystallites with c axes parallel to the applied field Ho 
(Ho 24 kOe), is more strongly broadened than the 
high-field side. Fits to the low-field region of the spec- 
trum yielded a crude estimate of the extra broadening, 
which becomes large at low temperatures as predicted by 
the disorder-driven NFL mechanisms discussed above. 

The magnetic susceptibility x(T) = M(H,T)/H, 
where M(H, T) is the bulk magnetization of the system, 
is strongly anisotropic in the Ce(Rui_ K Rh a ;)2Si2 series, 
with the c-axis susceptibility Xc(T) (Ho || c) much-larger 
than the a6-plane susceptibility Xab{T) (Ho _L c).E3 This 
suggests that the extra broadening observed for H || c 
might be due to disorder in the susceptibility similar to 
that found in UCu5_ x Pd x , and we were motivated pa 
measure the linewidth in a field-aligned powder sample.tj 
The powder was mixed with epoxy, which was allowed to 
harden in a magnetic field of 60 kOe. The torque on 



the anisotropic moment aligned the c axis of each single- 
crystal powder grain in the direction of the applied field 
before the epoxy hardened. 

The anisotropic susceptibility measured in this field- 
aligned powder sample is shown in Fig. ||. These data 
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FIG. 3. Temperature dependence of the anisotropic suscep- 
tibility in a field-aligned powder sample of CeRhRuSi2- Cir- 
cles: basal-plane susceptibility Xab(T). Triangles: c-axis sus- 
ceptibility Xc (T)- Solid curve: fit of single- ion Kondo-disorder 
model (Refs. jig and 0) to %c{T). Dashed curve: fit of Grif- 
fiths-phase model (Ref. [J) to Xc(T). 

agree well with measurements on a small single crystal of 
CeRhRuSi2 (not shown). A strong Curie- Weiss- like tem- 
perature dependence is found for Xc{T), whereas &iib{T) 
is small and only weakly temperature dependent .c3 The 
curves give fits of Xc(T) to the Kondo-disorder and 



Griffiths-phase models as discussed below in Sects. IV A 
and IV B| , respectively. 

Figure || also gives the 29 Si NMR spectra measured in 
the field-aligned powder sample for Ho || c and Ho _L 
c. It can be seen that, as expected, the line is wider 
for Ho || c than for Ho _L c. The small shoulder on 
the high-field side of the Ho || c line indicates that the 
alignment of crystallites in this sample is not perfect. The 
large linewidth anisotropy implies, however, that neither 
a small misalignment of the crystallites in the sample nor 
a slight misalignment of the sample with respect to Ho 
affect linewidth measurements appreciably for Ho || c. 
29 Si NMR spectra from a more completely aligned but 
smaller sample (not shown) confirmed this expectation. 
The misalignment does, however, preclude any attempt 
to obtain information about the shape of P(A) from the 
shape of the NMR line. 



III. SUSCEPTIBILITY INHOMOGENEITY AND 
NMR LINEWIDTH 

Since the NMR frequency shift of a given nucleus is 
determined by the interaction between its magnetic mo- 
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ment and those of the surrounding electrons, any spatial 
variation of the electronic magnetic susceptibility will be 
reflected in the NMR linewidth as a distribution of fre- 
quency shifts. A quantitative understanding of the sus- 
ceptibility inhomogeneity requires analysis of the relation 
between it and the NMR linewidth, independent of the 
particular mechanism which causes the inhomogeneity. 

The NMR frequency shift K measures the time- 
averaged effective field produced by the local moment 
at the resonating nucleus. In a paramagnet the relative 
shift Ki of the i th nucleus is related to the local suscepti- so that 
bility Xj associated with the j th /-ion electronic moment 

byka 



1/2 



Xj 



(1) 



where a,j is the hyperfine coupling constant between the 
j th moment and the i th nucleus. It is straightforward to 
carry out the spatial averages and show that 



K 



ax 



where a bar designates a spatial average in this and 
t he following. Similarly, the rms spread of shifts SK = 
(K 2 — X 2 ) 1 / 2 is related to the corresponding rms spread 
of susceptibilities Sx by 

5K = a*5 X , 

where a* is an effective hyperfine coupling constant, dis- 
cussed in more detail below. As a consequence 



6 X /x = 6K/(a*x)- 



(2) 



If each nucleus is coupled to more than one moment 
[cf. Eq. (|l|)], any spatial correlation between the mo- 
ment susceptibilities will affect the value of a*. There 
are two extreme limits in considering this spatial cor- 
relation. The term "long-range correlation" (LRC) will 
be used to describe the situation where the correlation 
length between local moments is much longer than the 
local-moment near-neighbor spacing. Similarly, "short- 
range correlation" (SRC) describes the situation where 
the variation of susceptibility from site to site is random 
or nearly so, i.e., where the correlation length which de- 
scribes this variation is of the order of or shorter than 
a lattice constant. Note that this correlation is only a 
phenomenological description of the inhomogeneous sus- 
ceptibility, and is not necessarily related to critical be- 
havior of the system. For a given system we do not know 
a priori which (if either) of these limits is applicable, al- 
though in the single-ion Kondo-disorder model we might 
expect that random ligand disorder would lead to rela- 
tively short-rangej-cprrelation. 

It can be shownllil that the values of a* in the LRC and 
SRC limits (denoted by and a| RC , respectively) are 
given by 



LRC 



^SRC 



X4 



Assuming for simplicity that the hyperfine coupling is 
predominantly to an effective number n e ff of /-ion near 
neighbors and is the same effective value a c g- for each of 
these neighbors, it follows that 



l SRC 



X SRC 



and 



l LRCI 



J LRC = n cffa c ff 



(3) 



In the LRC limit the fractional susceptibility inho- 
mogeneity Sx/x is given by the relative NMR line- 
width SK/ \K\. Since SK = a/Ho, where a is the rms 
linewidth in magnetic field units, we have from Eq. (|^) 

Sx _ SK __ SK 



X 



a LR C X 

a 



K H, 



(LRC limit) 



(4) 







Thus ct/(\K j Ho), which can be derived from the NMR 
data, is an estimator for Sx/x m the LRC limit. The 
corresponding estimator in the SRC limit can be obtained 
from <7/(|-?r| Ho) simply by scaling by the factor y/n^ 
[Eq- ©]: 

SX 



(5) 




(SRC limit) . 



The above assumes that the coupling constants a^- are 
not disordered, i.e., that they have the same values for 
crystallographically equivalent positions of nucleus i and 
/-ion j. If this is not the case and the are also disor- 
dered, then it can be showrH that 



SK 

W\ 



A 



1/2 



(6) 



where A is a term which expresses the effect of tb 
ordered a%j ■ Now in existing disorder-driven model: 
it is found that Sx/x varies considerably with % (with 
temperature an implicit parameter), tending to a value 
<; 1 at low temperatures and vanishing as x ~ y (high 
temperatures). On the other hand A is found to be in- 
dependent of x- Disorder in the a^- will therefore result 
in a nonzero intercept in a plot of SK/ \K\ vs. x, and 
its effect can be removed by subtracting this intercept in 
quadrature from the raw SK/ \K\ data. It can be shown 
that this correction is valid in both the LRC limit and 
the SRC limit. 

It should be stressed that this "NMR technology" is 
quite independent of the specific mechanism which causes 
the susceptibility inhomogeneity. 
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IV. DISORDER-DRIVEN NFL MECHANISMS 

In considering systems where NFL behavior is driven 
by disorder it is convenient to study the spatially- 
distributed local susceptibility Xj ■ Simple linear response 
theory shows that the zero-temperature local susceptibil- 
ity can be associated with a characteristic local energy 
scale Aj by 



1 

Xi <* XT 



(7) 



where Aj is essentially the excitation energy from the 
ground state to the first excited state. At finite tempera- 
tures T the local susceptibility x(A, T) depends strongly 
on the microscopic details which couple the magnetic de- 
grees of freedom. 

We can therefore speak of distributions of suscepti- 
bilities or energy scales, characterized by distribution 
functions P(A) and P(x), respectively; thus P(x,T) = 
P{A)/\d X (A,T)/dA\. Once P(A) [or P(x,T)} is known 
we can obtain spatial averages of physical quantities such 
as the n th moment x n (T) of the local susceptibility dis- 
tribution, which is given by 

POO 

X ir (T)= / X "(A,T)P(A)dA (8a) 



view of Eq. (0) one sees that regions of the system where 
Tk is very small (sites with very large low-temperature 
susceptibility) dominate the thermal and transport prop- 
erties. Miranda et aln have treated this picture in de- 
tail, and have shown that it predicts the observed low- 
temperature behavior of the Sommerfeld coefficient j(T), 
susceptibility x(T), and resistivity Ap(T) = p(T) — p(0) 
(7 oc x cx — lnT and Ap oc T, respectively) provided 
that P(T K ~>0) is finite. 

The resulting distribution function P(Tk) is given by 



P(T K ) = P(g) 



dg 



dT K 



g 2 P(g) 



(9) 



with g = 1/ \ti(Ef /Tk)- As a convenient parameteriza- 
tion of the Kondo physics we take the susceptibility to 
have the Curie- Weiss form 



x(T,T K )=C/(T + aT K ) 



(10) 



where C is the Curie constant. The value of a was es- 
timated by comparing-tliis Curie- Weiss law to the exact 
Bethe-ansatz solution;EZl the two functional forms differ 
by <10% for a » 2.5. Assuming a Gaussian distribution 
for P(g), the mean g and rms width Sg of the distribution 
can be found by fitting Eq. (j|) with n = 1 to the-Jiiea- 
sured bulk (i.e., spatially averaged) susceptibilityEJO 



X n P(x,T)d X . 



(8b) 



Knowledge of the first and second moments of P(x, T) is 
sufficient for interpretation of the bulk susceptibility and 
NMR linewidth data. 



A. Single-ion Kondo disorder 

We take the Ce-ion spins to be coupled to the 
conduction-electron spins by a s-f exchange coupling 
described by a coupling constant g = N(Ef)J, where 
N(Ep) is the density of conduction-electron states at the 
Fermi surface and J is the Ce-ion-conduction-electron 
exchange interaction. If the system is disordered on the 
ligand sites, as in CeRhRuSi2, g will be randomly dis- 
tributed according to a distribution function P(g). In 
the simplest picture of the Kondo effect, the Kondo tem- 
perature Tk , which characterizes the energy scale of the 
single-ion Kondo effect, is given by Tk — Efe^ 1 ^ 9 , 
where Ep is the Fermi energy. Thus a narrow distri- 
bution of g can lead to a wide distribution of Tk when g 
is small. In this picture we immediately identify A with 
T K - 

If the distribution function P(A) = P(Tr-) is broad 
enough so that P(Tk^0) does not vanish, then at any 
nonzero temperature T those / ions for which Tk < T 
are not compensated (i.e., are not described by Fermi- 
liquid theory) and give rise to the NFL behavior. In 



B. Griffiths-phase model 

In the Griffiths-phase model of NFL behavior! the 
low-energy physics is dominated by rare and large clus- 
ters which can tunnel over classically forbidden regions. 
These correlated regions are generated by above-average 
values of the RKKY interaction. The tunneling is pro- 
duced_by the spin-flip processes present in the Kondo 
effect In this scenario the Griffiths singularities ap- 
pear close to a QCP below percolation threshold and are 
therefore intrinsically related to QCP physics. It is intu- 
itively clear that the clusters can be effectively described 
in terms of two level systems, with tunneling energy E 
which is distributed over the sample due to the structural 
disorder. Obviously we have A = E in this picture. 

The distribution of E is obtained by mapping the prob- 
lem onto the Ising model in a transverse field; this proce- 
dure is valid in the limit of large magnetic anisotropy as 
appears to be the case in CeRuRhSi2 (cf . Fig. |J) . Then 
it can be shownEj that 




-l + A 



, < E < e ■ 
E > e , 



(11) 



where A is an exponent that determines the behavior of 
the response functions (0 < A < 1), and eo is a high 
energy cut-off which must be determined for each specific 
system. 
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As discussed above the local zero-temperature suscep- 
tibility in this picture is x(0, E) = C/E; large clusters 
with small energy scales have large susceptibilities. At 
high temperatures one expects the clusters to be disor- 
dered and behave paramagnetically, resulting in a Curie 
behavior x{E,T) = C/T for T > E. We therefore as- 
sume a Curie- Weiss interpolation formula 



X(E,T) 



C 



T + E 



(12) 



as for the Kondo-disorder model. But in the present 
case this approximation is intended to incorporate all the 
interaction effects which determine the susceptibility of 
a multi-ion cluster, not just single-ion Kondo physics. 

Using Eqs. (||), (|TT|), and ( [l2|) it is straightforward to 
sho-wEa that for T < e 



X(T) 



it AC /eoN 1 



and 



SX(T) 
X(T) 



eo 



1 — A 
2ttA 



m 



1/2 



A/2 



(13) 



(14) 



The critical behavior is determined by the single nonuni- 
versal temperature-independent exponent A. For A < 1 
the susceptibility diverges algebraically as T — > (H = 
0), and NFL behavior is obtained. The divergence in- 
creases (i.e., the NFL behavior becomes stronger) with 
decreasing A. The case A = 1 is marginal, and leads 
to logarithmic singularities as in the Kondo-disorder 
approach. 



C. Kondo disorder versus Griffiths singularities 

Both the Kondo-disorder and Griffiths-singularities 
pictures deal with a similar aspect of disorder, viz., the 
physics of rare events with large susceptibilities. It is 
clear, however, that the microscopic aspects of the two 
models are very different. The Kondo-disorder model 
uses non-interacting single- ion physics, and no aspect of 
the RKKY interaction is present. The Griffiths-phase 
approach, on the other hand, tries to take both RKKY 
and Kondo phenomena into account on an equal footing, 
and has a strong connection with QCP physics. 

From the point of view of local properties as measured 
in NMR spectra these approaches give similar results. 
The spatial properties of the two approaches are very 
different, however, since the formation of clusters in the 
Griffiths phase requires spatially extended structure. In 
this case one could look for the existence of clusters via 
superparamagnetic response, which is well understood in 
the context of spin glasses,E2l or for a momentum depen- 
dence of the inelastic neutron scattering. In addition, 
one would expect cluster formation to slow down the 
spin fluctuations relative to the free-ion fluctuation rate, 



which is essentially Tk- Experiments that are sensitive to 
the fluctuation rate may therefore be able to distinguish 
between the two theories. 



V. RESULTS AND DISCUSSION 



A. Bulk magnetic susceptibility 



A fit of the Kondo-disorder model result for x(T) 
[Eq. (|) with n = 1 and A = T K , and Eqs. (§) and 
@ for P{T K ) and x(T,T K )] to the experimental c- 
axis susceptibility Xc(T) is shown as the solid curve in 
Fig. ||. We obtain the same coupling constant distri- 
bution width Sg — 0.021 as Graf et al.$3 and a some- 
what smaller mean ~g = 0.160 compared to 0.175 from 
Ref. ^2|E3 The coupling constants are less widely dis- 
tributed—than in the NFL system UCus-^Pd^, x = 1.0 
and 1.5|1J consistent with weaker "nearly NFL" behavior 
in CeRhRuSi 2 . 

Figure |^ also gives the Griffiths-phase model predic- 
tion for Xc(T) (dashed curve), obtained by fitting y(T) 
from Eq. (||) (n = 1), with A = E and using Eqs. ( fil| ) 
and (|l2|), to the bulk c-axis susceptibility. The best fit is 
given by the dashed curve in Fig. |[ It can be seen that 
at low temperatures the Griffiths-phase fit curve overes- 
timates the experimental data slightly. This is to be ex- 
pected, since in the simple Griffiths-phase model there is 
no possibility of a return to Fermi-liquid behavior at low 
temperatures: the system is a true NFL as long as A < 1. 
But the susceptibility data begin to exhibit the satura- 
tion expected from the conclusions of Graf et aZ.,t3 and 
therefore are not well described by an algebraically diver- 
gent temperature dependence. There is, however, a re- 
gion of intermediate temperatures in which both Kondo- 
disorder and Griffiths-singularity models agree very well 
with experiment. From the Griffiths-phase fit in this re- 
gion we obtain e = 170 ± 10 K and A = 0.88 ± 0.02. The 
latter value is considerably larger (i.e.^-jthe NFL behavior 
is weaker) than found in UCus-^PdJj as in the Kondo 
disorder model. 

Figure || shows the distribution functions P(A) which 
result from the Kondo-disorder (A = Tk) and Griffiths- 
phase (A = E) model fits. It can be seen that the two 
functions are very different. The Kondo-disorder distri- 
bution function P(Tk) exhibits a maximum near 12 K 
and is small below ~1 K and above ^100 K, whereas 
the Griffiths-phase distribution function P(E) is broader 
and diverges weakly as E — > 0. These differences are 
much more marked in P(A) than in the corresponding 
fits to the susceptibility (Fig. ||); spatially-averaged ex- 
perimental quantities are insensitive to the exact form of 
P (A). It is clear that the Griffiths-phase model could fit 
the data better if P(E) — » as E — ► 0, which in an ex- 
tended Griffiths-phase picture would occur if there were 
an upper cutoff on the cluster susceptibility. 
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FIG. 4. Distribution functions -P(A) of characteristic en- 
ergies A in CeRhRuSi2, obtained from fitting the Kon- 
do-disorder (A = Tk) and Griffiths-phase (A = E) 
models to the bulk c-axis susceptibility (Fig. ||). Solid 
curve: Kondo-disorder distribution function P(Tk) [Eq. ([])]. 
Dashed curve: Griffiths-phase distribution function P(E) 
[Eq. (El)]. 



Thus CeRhRuSi 2 does not exhibit "true" NFL be- 
havior. In the Kondo-disorder model, which shows this 
most explicitly, the best fit indicates that all spins are in 
a Kondo-compensated Fermi-liquid state at low enough 
temperatures. This is consistent with the results of Graf 
et alH that j(T) -> const, and Ap(T) cx T 2 as T 0. 
We note again, however, that as mentioned in Sec. | the 
saturation of 7(T) does not necessarily indicate Fermi- 
liquid behavior at low-itemperatures; other physics, such 
as magnetic freezing j£d'E3 may be involved. 



B. 29 Si NMR linewidths 

The 29 Si c-axis NMR shift K c and linewidth a c are 
plotted against the c-axis bulk susceptibility Xc in Fig. ||, 
with temperature an implicit parameter. (The a6-plane 
parameters K a b and a a b, not shown, are small and only 
weakly temperature dependent.) It can be seen that a c 
varies more rapidly with Xc than K c , as expected quality 
tively from disorder-driven theories of NFL behavior .ETt2l 
Although the expected linear relation between K c and 
Xc is observed at high temperatures (small X c)> Kc{Xc) 
tends to a constant for large Xc- This saturation is not 
well understood, but may be due to a small amount of 
second phase; this could have a strong Curie- Weiss-like 
susceptibility but little effect on the NMR shift since the. 
number of nuclei in the second phase would be smalLEil 
The observed nonlinearity is not more than ~20% and 
does not affect our conclusions significantly. 

Figure | plots a c /(K c H ) = SK C /K C = 6K c /(a£ KcX c) 
[Eq. (0)] versus X c- As discussed in Sec. Ill a c /(K c H ) 
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FIG. 5. 29 Si c-axis NMR shift K c (circles) and linewidth c c 
(triangles) vs. c-axis susceptibility Xc in CeRhRuSi2, with 
temperature an implicit parameter. Spectrometer frequency 
20.220 MHz (applied field ~ 23.9 kOe). The linewidth a c 
varies more rapidly with \c than K c , as expected from disor- 
der-driven theories of NFL behavior (Refs. §-|lo|). 



is an estimator for 5x/x in the LRC 
extrapolate to a non-zero value as Xc 



limit. The data 
— * 0, which indi- 



cates that the coupling constant is also disordered. 
We therefore subtracted the extrapolated intercept from 
the raw values in quadrature [cf. Eq. (|g)] to obtain cor- 
rected data in the LRC limit, also shown in Fig. || (tri- 
angles). These corrected data represent 5K C / (a^ RC Xc) 
with SK C due only to susceptibility inhomogeneity. The 
corresponding (corrected) values of SK C / '(igRcXc) were 
obtained from Eq. (||) of Sec. Ill, with n c s chosen as 
described below. 

For both the Kondo-disorder and Griffiths-phase mod- 
els Sx/x was calculated from Eq. (||) (n — 2) and x(T) 
with no further adjustable parameters, since P(A) had 
been previously determined by the fits to the bulk sus- 
ceptibility. Figure |?] compares 6K C / (a*Xc) in both limits 
with the theoretical behavior of Sx/x from these theories, 
again with temperature an implicit parameter. It can be 
seen that the theoretical predictions are similar and that 
they both overestimate SK C / (a^ c Xc) considerably, but 
that the agreement with SK c /(a^- RC Xc) is excellent when 
n c ff in Eq. (||) is taken to be 6. 

That this value is sensible can be concluded from 
examination of the A^Ba-type crystal structure of 
CeRhRuSi2, shown in Fig. |[ where it can be seen that 
each Si site is coordinated by four Ce nearest neighbors 
and one Ce next-nearest neighbor. Thus n c ff is approx- 
imately the coordination number for the first two near- 
neighbor shells and is therefore reasonable, given the ap- 
proximation of an effective number of equally-coupled 
neighbors. 

For CeRhRuSi2 we do not have the independent veri- 
fication of the SRC limit that was available from com- 
parison of NMR and muon spift, rotation (/iSR) spec- 
tra in the case of UCu5_ 2; Pd x .E 2 J (For a review of the 



7 



0.4 



0.3 



CeRhRuSL 



^ 0.2 



0.1 
0.0 



T • 

• fx 



a 

corrected 



20 



40 



60 



Susceptibility (10 emu/mole) 



1.0 



« 0.5 

CO 



0.0 



CeRhRuSi 2 


ii i i 


Kondo Disorder 




Griffiths Phase 




T> 


SRC limit 








• ■ 

LRC limit 



20 



40 



60 



Susceptibility (10 emu/mole) 



FIG. 6. Plot of relative spread a c /(K c Ho) of c-axis 29 Si 
NMR shifts in LRC limit versus c-axis susceptibility \ c in 
CeRhRuSi2- See text for symbol definitions. Circles: raw 
data before correction for disorder in the hyperfine coupling 
constant. Triangles: data corrected for coupling constant dis- 
order (see text). 



/j,SR technique see, for example, Ref. [53].) The latter 
alloys have a cubic crystal structure and the positive- 
muon (/U + ) interstitial stopping sites possess octahedral 
and tetrahedral point symmetries, which are sufficiently 
high to render the fi + frequency shift isotropic. Then 
the /i + linewidth reflects the susceptibility inhomogene- 
ity rather than anisotropic powder-pattern broadening. 
Preliminary /iSR measurements in an unaligned pow- 
der sample of CeRhRuSi2u show that in this alloy the 
anisotropic contribution dominates the powder-pattern 
linewidth, much as in the unaligned-powder spectrum of 
Fig. |^, and the disorder-induced broadening cannot be 
determined accurately. 

Unfortunately field-aligned powder samples cannot be 
used in [iSR experiments. The packing fraction of the 
powder must be small 20%) in order to allow free ro- 
tation of the powder grains during alignment, and then 
only a correspondingly small fraction of the muons stop in 
the sample; the rest stop in the epoxy and give a spurious 
signal. Thus we cannot confirm the SRC limit by com- 
paring results between NMR and /iSR. We also note that 
no other nucleus in CeRhRuSi2 is favorable for NMR; sta- 
ble Ce isotopes possess no nuclear magnetic moment, and 
Ru and Rh isotopes have very small gyromagnetic ratios. 
Nevertheless, the SRC-limit estimate of $Xc/Xc is in ex- 
cellent self-consistent agreement with the disorder-driven 
models. 



VI. CONCLUSIONS 

The picture that emerges from our 29 Si NMR stu- 
dy of CcRhRuSi2 exhibits similarities and differences 
when compared to the preceding NMR investigation of 



FIG. 7. Experimental and theoretical relations between 
rms susceptibility inhomogeneity and bulk susceptibility in 
CeRhRuSi2. Data points: c-axis NMR estimator 5K c /(a Xc) 
of relative rms susceptibility spread &x/X- Circles: LRC limit 
(see text). Triangles: SRC limit. Solid curve: 5x1 X i n the 
single-ion Kondo-disorder model. Dashed curve: &xlx m the 
Griffiths-phase model. Agreement is found between predic- 
tions of disorder-driven NFL theories and data in the SRC 
limit for effective nearest-neighbor number n e s = 6. 



UCus-zPda;, x = 1.0 and 1.5. BO The most important 
similarity is the fact that in both cases the NMR data 
are in excellent agreement with predictions of disorder- 
driven theories of NFL behavior. Gap results therefore 
confirm the conclusions of Graf et alS3 that such a mech- 
anism drives NFL properties in CeRhRuSi2- The most 
important differences between the two systems are that 
in CeRhRuSi2 (1) within the single-ion Kondo-disorder 
model the disorder is not enough to prevent a return to 
Fermi-liquid behavior at temperatures < 1 K, and (2) the 
determination of the appropriate correlation length limit 
(LRC or SRC) has not been made independently of com- 
parison with theory. The agreement between theory and 
experiment assuming the SRC limit (Fig. is satisfac- 
tory. 

From the experimental point of view, the relatively 
small differences between the predictions of the single- 
ion Kondo-disorder picture and the Griffiths-phase the- 
ory show how difficult it is to discriminate between these 
two mechanisms for disorder-driven NFL behavior in the 
NMR linewidths. We speculate, however, that the dy- 
namics of the spins will be quite different in the two cases, 
particularly at low temperatures. 

The single-ion Kondo disorder model predicts inhomo- 
geneous relaxation due to the distributed Tk- This mech- 
anism yields a spatially averaged spin-lattice relaxation 
rate T^ X (T) = J dTjc P (Tk ) 1 (T, Tk ) ■ It is convenient 
to approximate T^~ x (T,Tk) by 



[ 1/T K , T>T K 
Tf (T, T K ) oc < ' ' 
1 T/T* , 



T < T, 



K 



8 



Rh, Ru 




Ce( 1 



FIG. 8. Crystal structure of Ce^hi-^Ru^^. Each Si 
site is coordinated by four Ce nearest neighbors (no. 1) in the 
adjacent basal plane and one Ce next-nearest neighbor (no. 2) 
in the opposite basal plane. 



which captures the crossover to Fermi-liquid (Korringa) 
behavior for T < Tk- For a model rectangular P(Tk) 
given by 



1 



P(T K ) = { T M - T„ 
0, 



T m < Tk < Tm ■ 
otherwise, 



where T m and Tm are minimum and maximum values of 
Tk, respectively, it is straightforward to show that Tj -1 
varies linearly with T for T < T m and goes smoothly to a 
constant for T > Tm- Such behavior would be generally 
expected to characterize T^ 1 as long as P(Tk— >0) = 0. 

Thus in this scenario Tj~ depends on temperature but 
is independent of resonance frequency lu. 

In contrast, it can be easily shown from the expression 
for the dissipative dynamic susceptibility x"(w) in the 
Griffiths-phase theorya 

X"{u) oc uj- 1+x ta,nh(huj/k B T) 

that, assuming the validity of this picture for the very 
low nuclear (muon) frequencies (Hlu <C ksT) 



TV 1 



-l + A 



independent of temperature. The frequency is given by 
uj = "/Ho, where 7 is the nuclear (muon) gyromagnetic 

moment. Thus the dependence of T^ 1 on temperature 
and Hq differs consierably between the two theories. 

It should be noted that the NMR linewidths obtained 
from experimental data are not necessarily true sample 
rms averages, if the lines have extended shoulders which 
are lost in the noise and not taken into account. The fact 
that good fits are obtained with Gaussian lines seems to 
make this unlikely, but it is not hard to see that at low 



temperatures P(x) should be broad and asymmetric in 
both the Kondo-disorder and Griffiths-phase models. If 
the experimental linewidth characterizes only "typical" 
environments it will underestimate the true sample av- 
erage. This would render our quantitative results some- 
what uncertain, but would not invalidate the conclusion 
that disorder is an important element in the NFL behav- 
ior of CeRhRuSi2 above 1 K. 

Finally, we discuss the relation of the disorder- 
driven theories to the observed crossover to a new 
regime (Fermi-liquid behavior, cluster formation, mag- 
netic freezing, etc.) in CeRuRhSi2 below 1 K. The 
crossover is described empirically by the Kondo-disorder 
model, which by itself gives no clue as to why there should 
(or should not) be a suppression of low Kondo temper^, 
atures. Recently, however, Miranda and DobrosavljvicS 
have reported a microscopic calculation of the form of 
P(Tk) for various levels of disorder. They find that the 
distribution is singular only for sufficiently strong disor- 
der, whereas for slightly weaker disorder P(Tk) — * at 
small Tk, onsistent with a return to Fermi- liquid behav- 
ior at the lowest temperatures. This feature is in at least 
qualitative agreement with our results. 

To explain the crossover the Griffiths-phase scenario 
would have to be extended beyond its simplest form 
to include a mechanism which reduces the response of 
the largest clusters. The mechanism behind such a re- 
duction could be the breakdown of the assumption of 
strong single-ion anisotropy made in the Griffiths-phase 
theory;Ej transverse fluctuations of the Ce ions might 
constitute a damping mechanism which rounds off the 
Griffiths singularities. Alternatively, superparamagnetic 
spin freezing . o£ t he clusters could occur at very low 
temperatures .EdE3 

It is possible that similar crossovers occur in other NFL 
materials, perhaps at temperatures which have not yet 
been explored. In any event, our experimental findings 
indicate that further development of current theories of 
disorder-driven NFL behavior is required to understand 
NFL phenomena at low temperatures. 
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